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ABSTRACT 

X-ray images of galaxy clusters and gas-rich elliptical galaxies show a wealth of small-scale features 
which reflect fluctuations in density and/or temperature of the intra-cluster medium. In this paper we 
study these fluctuations in M87/Virgo, to establish whether sound waves/shocks, bubbles or uplifted 
cold gas dominate the structure. We exploit the strong dependence of the emissivity on density 
and temperature in different energy bands to distinguish between these processes. Using simulations 
we demonstrate that our analysis recovers the leading type of fluctuation even in the presence of 
projection effects and temperature gradients. We confirm the isobaric nature of cool filaments of gas 
entrained by buoyantly rising bubbles, extending to 7' to the east and south-west, and the adiabatic 
nature of the weak shocks at 40" and 3' from the centre. For features of ~ 5 — 10 kpc, we show that 
the central 4' x 4' region is dominated by cool structures in pressure equilibrium with the ambient 
hotter gas while up to 30% of the variance in this region can be ascribed to adiabatic fluctuations. 
The remaining part of the central 14' x 14' region, excluding the arms and shocks described above, is 
dominated by apparently isothermal fluctuations (bubbles) with a possible admixture (at the level of 
~30%) of adiabatic (sound waves) and by isobaric structures. Larger features, of about 30 kpc, show 
a stronger contribution from isobaric fluctuations. The results broadly agree with an AGN feedback 
model mediated by bubbles of relativistic plasma. 


1. INTRODUCTION 


It is well known that the cores of relaxed clus¬ 
ters are strongly perturbed by the outflows of rel¬ 
ativistic plas ma from a cen tr al su p ermassive black 


2000; McNamara et all 200(1 Fabian et all 

2000; 

David et all 200i[|Nulsen et all 20021; Birzan et al.l 

2004; 

Dunn & Fabian 2000), suggesting that the amount of 


energy supplied by the AGN is sufficient to offset gas 
cooling losses. In the vicinity of the AGN these outflows 
inflate “cavities/bubbles” in X-ray-emitting gas, which 
evolve under the ac t ion of a buoyancy force (e.g., 
IChurazov et al.l 12001b iRevnolds et all 120051) . Among 
plausible channels which transfer energy from the out¬ 
flow to the g as two processes rec e ived much attention : 


shocks 


(e.g, 


hrbi an et all 12003)1 iForman et all 120071 


I Randall et al.l 120111) and purely subsonic interaction 


of the rising bubbles through the X-ray atmosphere, 
which uplifts cool gas, drives turbulen ce and excites 
gravity waves (e.g. IChurazov et al.l 120021) . When these 
processes are operating at the same time problems arise 
in differentiating among them in order to assess their 
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relative roles. In this paper we present a method 
to broadly identify the type of process, either adia¬ 
batic, isobaric or isothermal, that leads to the largest 
amplitude fluctuations. 

Here we use a simple technique (C hurazov et al. in 
prep., see also iZhuravleva et all I2015T) of combining in¬ 
formation contained in X-ray images in different energy 
bands to determine which process makes the largest con¬ 
tribution to the observed fluctuations of the intracluster 
medium (ICM) in M87/Virgo. We do this by calculating 
the cross power spectrum for two Chandra images and 
calculating the correlation coefficient and relative ampli¬ 
tudes of fluctuations in several energy bands, namely 0.5- 
1 keV, 1-3.5 keV and 3.5-7.5 keV. M87/Virgo is a canoni¬ 
cal nearby example of a cool core cluster with prominent 


substructure (e.g-, Bohringer et al.l 

1993; Forman et al. 

2005, 2007j; Simionescu et akl 2010; 

Million et al. 201(3). 


Throughout the paper we adopt a distance to M87 of 16 
Mpc, V corresponds to 4.65 kpc. Its proximity and rel¬ 
atively low gas temperature makes it an ideal object for 
study based on Chandra X-ray images. In an accom¬ 
panying paper we also apply this analysis to the Perseus 
cluster (Zhuravleva et al. in prep). 

In Sec. [2] we detail the effect that density fluctuations 
have on X-ray images in different bands and how this can 
be exploited to differentiate between isobaric, adiabatic 
and isothermal processes. In Sec. [3] we introduce the 
power and cross spectrum measures that will be used to 
characterize the images and we test these statistics with 
simulated cluster images in Sec. [4] We proceed in Sec. [5] 
with real cluster data, using Chandra images of M87 to 
determine the dominant type of fluctuation in different 
regions of the ICM. We interpret these results in term 
of the energy budget of the fluctuations in Sec. [HI and 
summarize our findings in Sec. [7] 
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2. ENERGY DEPENDENT IMPACT OF FLUCTUATIONS 
ON X-RAY IMAGES. 

Many well known processes can cause fluctuations of 
the ICM density and temperature, which in turn give 
rise to perturbed X-ray images. Among them: (i) 
weak shocks and sound waves, (ii) bubbles of relativistic 
plasma that make X-ray cavities in the gas, (iii) gas dis¬ 
placement from its equilibriumposition by buoyant bub¬ 
bles and subsonic gas motionifj Each of these processes 
will lead to specific energy dependent signatures in the 
X-ray data. For instance, weak shocks and sound waves 
cause density fluctuations that are positively correlated 
with temperature fluctuations, as long as thermal con¬ 
duction is suppressed. On the other hand, cavities asso¬ 
ciated with the bubbles of relativistic plasma, may have 
similar depressions in X-ray images in all bands (if the X- 
ray emissivity of the plasma inside the cavities is small). 
Finally, subsonically displaced gas lumps may conserve 
their entropy (once again, provided that conduction is 
suppressed) and remain in pressure equilibrium with the 
ambient gas. To maintain pressure balance, these lumps 
will show anti-correlated temperature and density fluc¬ 
tuations. Thus, even if the true ICM equation of state 
is adiabatic with 7 = 5/3, the observed correlations of 
the density and temperature fluctuations may not cor¬ 
respond to the 7 = 5/3 equation of state. Below we 
classify different types of fluctuations based on their ap¬ 
parent equation of state, based on the relative amplitudes 
of fluctuations present in X-ray images in different energy 
bands. 

As mentioned above, the three aforementioned pro¬ 
cesses differ by the temperature changes associated with 
a given amplitude of the density fluctuation. Our goal is 
to use the images in different energy bands to determine 
the type of process which gives rise to the density fluctu¬ 
ations. Assuming that the flux emitted per unit volume 
is F = p 2 Ab(T), where p is the density and Aj,(T) is 
the emissivity in band b as a function of temperature 
T, all that remains is to calculate the relation between 
fluctuations of p and T for each type of process and the 
emissivity Ab(T) for each energy band. 

We estimated Af,(T) by folding an XSpec APEC model 
through the response function of Chandra’s ACIS-I de¬ 
tector, and calculate the flux in each band as a func¬ 
tion o f tem perature, following the procedure described in 
iForman et al.1 (|2007t ). For this study we use two energy 
bands: 1.0-3.5 keV and 3.5-7.5 keV (“soft” and “hard” 
bands, respectively). The choice of the reference energy 
bands is driven by two factors. On the one hand, using a 
combination of images in very soft and very hard bands 
maximizes the sensitivity to the temperature variations. 
On the other hand, we need sufficient numbers of counts 
in each image to produce meaningful results. In Chura- 
zov et al. (in prep) we address this issue and show that 
the pair we use is close to optimal for the temperature 
characteristic of M87. 

For unperturbed temperature and density To and po, 
the perturbed values are 


7 This list is by no means complete, but for the sake of simplicity 
we concentrate on these three processes. 


1. T = To(p/po) 1 for an isobaric fluctuation, 

2. T = T 0 (p/po) 2 / 3 for adiabatic fluctuations and 


3. T = To for isothermal fluctuations. 

Therefore, the change of the flux 5Fb in band b in re¬ 
sponse to a change in density 


p = G x po 


( 1 ) 


is 


6 Ft 

~Fb 


b_ _ p 2 A b (T)-p 2 0 A b (T 0 ) _ p 2 0 G 2 A b (T0G’ 1 ) _ 

PoA b (To) PoA|,(Tb) 


G A b (T 0 G n ) 
ATfTbl 


- 1 


( 2 ) 

(3) 


where the exponent p is equal -1, 2/3 and 0 for iso¬ 
baric, adiabatic and isothermal processes, respectively. 
In the limit of small density fluctuations G = 1 + x, 

SFb 

where x <C 1 , the above expression reduces to —— = 

kb 

d log A b \ 

p ——— I x. 


dlogT J 

In FigO] we show the ratio of relative amplitudes 
Ty 1 / 7 ?= °f the fluctuations in the hard and soft band 
images corresponding to the same fluctuation in den¬ 
sity. Adiabatic and isobaric cases are shown in the top 
and bottom panels, respectively. Isothermal fluctuations 
have equal amplitudes in all bands, i.e. they produce a 
ratio equal to 1, irrespective of temperature. The green 
thick line in each set corresponds to the limit of small 
amplitude fluctuations. To illustrate the nonlinear ef¬ 
fects of large amplitude fluctuations we also show the 
curves for G = -F = 0.5 (red, solid line), G = 1.5 (blue, 

dashed line) and G = 2.0 (black, dot-dashed line). In 
the temperature range relevant to M87, 1.3-3 keV, the 
three types of process are clearly distinguishable, pro¬ 
ducing larger, smaller or equal amplitude fluctuations in 
the hard band than in the soft band for adiabatic, iso¬ 
baric and isothermal fluctuations, respectively. Factor of 
~2 fluctuations in density shows substantial deviations 
from the small amplitude limit and marked asymmetry 
between positive and negative fluctuations (G = 2 or 
G = 0.5). For smaller amplitudes (less than ~50%), the 
green curves provides a reasonable approximation of the 
flux ratio. 

One noticeable feature of isobaric fluctuations is that 
around 2 keV, which is close to the average temperature 
in M87, the ratio is close to zero. This is be¬ 

ef log A 


cause, for the hard (3.5-7.5 keV) band, 


T 


SF h 


dlogT 


2 at 


2 keV and so —— ~ 0. The 3.5-7.5 keV image is 
Fh 

therefore a good proxy for densit y flu ctuations produced 
by non-isobaric fluctuations (|Forman et al.ll2007 ). i.e., it 
is only sensitive to adiabatic and isothermal fluctuations. 
In contrast, for the 1.0-3.5 keV band the emissivity is 
essentially independent of temperature, as noted before, 
dlog A 


Thus 


dlogT 


0 and 


2x. Therefore this band is 
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Fig. 1.— Ratio of relative amplitudes of flux fluctuations in the 
“hard” (3.5-7.5 keV) and “soft” (1.0-3.5 keV) bands as a function 
of initial temperature To, for adiabatic processes (top panel) and 
isobaric processes (bottom panel). A factor G ( see EqfT) charac¬ 
terizes the amplitude of the density fluctuation. The thick green 
line in each plot corresponds to a limit small density fluctuations 
(G — 1 <C 1), while the thin lines correspond to fluctuations with a 
given finite amplitude, as labeled in the plot. For the typical gas 
temperatures in M87 (T ~ 2 keV) the amplitudes for isobaric and 
adiabatic processes differ strongly, opening the possibility to use 
this difference to distinguish these processes in Chandra images. 

a proxy for density variations independently of the type 
of process that produces them. Figure [2] shows a soft 
(0.5-3.5 keV) and a hard (3.5-7.5 keV) image. The arms 
disappear in the hard band image, consistent with an 
isobaric nature interpretation, while the rings appear as 
larger amplitude surface brightness fluctuations in the 
hard band than in the soft band image, as expected for 
adiabatic weak shocks. 

As is clear from Fig. |T] the ratio for isobaric 

fluctuations changes sign around 2 keV and is also sen¬ 
sitive to the amplitude of fluctuations. Therefore one 
might expect to find either correlation or anti-correlation 
of fluxes in these two bands, although the absolute value 
of this ratio for isobaric fluctuations remains significantly 
smaller than for isothermal or adiabatic fluctuations. 

3. POWER SPECTRA 

We now proceed with the analysis of scale-dependent 
fluctuations, including faint ones, which we would like to 
characterize in a statistical sense through the power spec¬ 
tra and cross power spectra of two Chandra images. We 
will focus on fluctuations around the large scale intensity 
profile of the cluster by first dividing the images by the 
projected spherically symmetric /3-model which best fits 
the azimuthally-averaged radial profile. The fluctuation 
images I = — 1 in each energy band are used in the 

analysis below. Here I ra w is the original image and Ip is 
the symmetric /3-model. The values of the slope /3 and 
core radius r c of the /3 models fitted to the image in each 


band are given in Sec. [5] 

We will use the method described in lArevalo et al.1 
(|2012T) to filter the fluctuation images in different spa¬ 
tial scales in order to calculate the power and cross spec¬ 
tra. This method corrects for gaps and irregular bound¬ 
aries in the images. The 2D power spectrum can be 
converted into a 3D power spectrum_through th e de - 
projection procedure explained in lChurazov et al.l (120121 ) 
and finally converted into the amplitude of fluctuations 
in 3D as = i JP^d * 47rfc 3 /2, where k is the inverse 
of the spatial scale in units of arcsec -1 , and P^d is the 
power per mode in the 3D power spectrum of emissivity 
fluctuations, reconstructed from the observed 2D power 
spectrum. Note that we use k = l/l where l is the spa¬ 
tial scale, without the usual factor of 2n in this relation. 
The factor 2 in the above expression is used to recast the 
amplitude in terms of linear density fluctuations instead 
of volume emissivity variations, since X-ray flux scales 
as oc p 2 and for small fluctuations this approximates to 
p 2 /p 2 = (1 + x ) 2 ~ 1 + 2x where p 0 is the unperturbed 
density and x is the density fluctuation. The cross spec¬ 
tra Pi ,2 are calculated using a similar procedure, by fil¬ 
tering the soft and hard band images with the Mexican 
hat filter, multiplying the filtered images in two bands by 
each other and calculating the mean value over all pixels. 

In the discussion below, we assume that fluctuations, 
associated with adiabatic, isobaric and isothermal pro¬ 
cesses are spatially uncorrelated (in 3D). This assump¬ 
tion simplifies the interpretation of the results by re¬ 
moving all cross-terms for different types of fluctuations. 
While one can easily imagine a situations when different 
types of processes are correlated (e.g., buoyant bubbles 
are often accompanied by cool entrained gas), the proce¬ 
dure described below should be able to establish which 
type of process dominates the fluctuations in the images 
in a statistical sense. 


3.1. Correlation coefficient 

Since the relative amplitudes of the soft and hard band 
fluctuations differ for different processes, fluctuations ob¬ 
served in two energy bands are expected to perfectly 
match each other (i.e. I\ = a/ 2 , where a is a constant) 
only if the observed fluctuations are dominated by one 
particular type of process. To this end, a handy quantity 
is the correlation coefficient, which verifies if fluctuations 
in one band are linearly related to the fluctuations in an¬ 
other band. Given two fluctuation images, the correla¬ 
tion coefficient Cm can be computed as 


h x h 


(4) 


where the subscript M refers to the measured value as 
opposed to the predicted value. In practice, we calculate 

a scale dependent correlation coefficient Cm = . 1,2 , 

VP1P2 

where Pi, P 2 and P i2 are the power spectra of / 1 , I 2 
and their cross spectrum, respectively. The contribution 
of the Poisson noise to Pi, P 2 is removed by generating 
several images of Poisson fluctuations around the fitted 
/3-model, calculating their power spectra and averaging 
the results. The averaged Poisson noise power spectrum 
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Fig. 2.— 0.5-3.5 keV and 3.5-7.5 keV images of M87 divided by their be st-fitting /3-models. The 0.5-3.5 keV image shows complicated 
substructure, while the 3.5-7.5 keV looks strikingly different. As shown by I Forman et ahl 1120071) , for gas temperatures, characteristic for 
M87, the surface brightness in this band is insensitive to isobaric fluctuations, but reflects projected pressure fluctuations. The lack of 
prominent “arms” in this image suggests that they are in pressure equilibrium with the ambient gas. 


is then subtracted from Pi and Pi- Unlike the standard 
definition of the correlation coefficient in signal process¬ 
ing, we keep the sign of correlation coefficient, i.e. the 
allowed range of Cm is between -1 and 1, to differentiate 
correlated from anti-correlated fluctuations. 

For a single type of process with small amplitude fluc¬ 
tuations in an isothermal cluster, the image in one band 
is just a linear transformation of the image in the other. 
Therefore, the correlation coefficient should be 1 or -1. In 
reality, clusters might be perturbed by different processes 
with comparable amplitudes. Some of these will have 
larger responses in the hard band than in the soft band 
and some other, spatially uncorrelated fluctuations, will 
have the opposite effect, so that the final images in both 
bands might also appear uncorrelated. This mixture of 
processes will be reflected in the measured correlation 
coefficient. 

For a combination of different (independent) processes, 
the fluctuation image in band b can be expressed as 
lb oc Y^, a b,iXi , where X, is proportional to the line-of- 
sight emission measure fluctuation relative to the un¬ 
perturbed model, associated with the process i, and 

Oib,i = ^2 + Vi d d°ogT ) the response in band b to this 

fluctuation. The expected correlation coefficient between 
soft and hard band images is then: 


C F = 


J2i °l,i° 2,i(x?) 


(5) 


where the cross terms ( XiXj) have been dropped by the 
assumption of uncorrelated variations. It is clear from 
this equation that if a single process dominates the den¬ 
sity fluctuations, e.g. (X?) » {Xj) for all j and all 
aq t i and i are of the same order, then the correlation 

coefficient converges to Ce = } i. e ., to a value 

of 1 or -1. On the other hand, for a mixture of processes 
with comparable amplitudes, the correlation coefficient 
takes a value \Ce\ < 1, that can be predicted once the 
values of a are known. 

Notice, that Cm will be poorly constrained if the lead¬ 
ing process does not produce a strong signal in one of 
the bands (i.e. if the response in I\ or I 2 to this process 
is close to zero). For T = 2 keV gas, this is the case for 
isobaric processes in the 3.5-7.5 keV band (see Fig. |IJ. 
As is clear from this figure, even for small variations of 
temperature around 2 keV, the sign of the ratio of ampli¬ 
tudes changes, leading to large uncertainty in the value 
of the correlation coefficient. 

3.2. Relative amplitudes of soft and hard band 
fluctuations 

Another useful quantity is the relative amplitudes of 
fluctuations in two energy bands R. The correlation coef¬ 
ficient C defined above in Eq. [^equals the cross spectrum 
normalized by the variance of both bands. To estimate 
the relative amplitude of fluctuations, R 1 we normalize 
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the cross spectrum by the variance of the soft band only, 

R= 1\ =CX ^' (6) 

R is proportional to the correlation coefficient and to the 
ratio of the amplitudes of the fluctuations in the hard 
and soft bands. The fact that we multiply the simple ra¬ 
tio of powers by the correlation coefficient preferentially 
highlights the ratio of correlated (or anti-correlated) fluc¬ 
tuations, which are the ones revealing we wish to study. 
We use the form -2^1 to calculate this quantity rather 

than C x ^5 because the first expression removes much 

of the uncertainty related to subtracting Poisson noise, 
which is high in the hard band. Since this noise is uncor¬ 
related in the two bands, the Poisson noise cancels out in 
the cross term Pi , 2 , and the denominator only includes 
the soft band power Pi which has a lower Poisson noise 
level. 


3.3. R — C map 

In the case of two or three types of independent fluctu¬ 
ation processes being comparably important, we expect 
that the relative amplitudes and the correlation coeffi¬ 
cient of the soft and hard band images will have inter¬ 
mediate values. 

Figure [3] shows the expected correlation coefficient 
(magenta contours) and amplitude ratio (blue contours) 
for the simulated images in the 1-3.5 keV and 3.5-7.5 keV 
bands. The top figure is calculated for an unperturbed 
temperature To = 2 keV and the bottom figure is for 
To = 1-6 keV. We normalize the total variance of density 
fluctuations, i.e., the sum of the variances correspond¬ 
ing to isobaric, adiabatic and isothermal processes, to 1. 
Therefore, in 3D Cartesian space, where the amplitude 
of each type of fluctuation corresponds to one of three or¬ 
thogonal axes, the parameter space considered here cor¬ 
responds to one octant of a spherical shell, which we 
project into the (x,y) plane. The amplitude of adiabatic 
fluctuations increases towards the right, isobaric fluctua¬ 
tions upwards and isothermal fluctuations decrease with 
increasing distance to the origin z = \/l — ( x 2 + y 2 ), so 
that pure adiabatic fluctuations are mapped into the bot¬ 
tom right corner, pure isobaric in the top left and pure 
isothermal to the origin. The correlation coefficient and 
amplitude ratio contours form an almost orthogonal pat¬ 
terns in these R— C maps. Therefore, for our idealised 
model of three uncorrelated processes, it is possible to 
find the coordinates of the intersection of the relevant R 
and C contours and from these to determine the fraction 
of the total variance produced by each type of density 
fluctuation. These unperturbed temperatures, Tq = 1.6 
and To = 2 keV were chosen to bracket the average tem¬ 
perature in M87 within 6' from the centre, which is the 
area that will be studied below. For the real cluster, or 
simulations calculated using the real cluster temperature 
radial profile, the measured R and C coordinates should 
lie between the limiting maps shown in Fig. [3] 
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Fig. 3.— R-C maps: Predicted correlation coefficient C (magenta 
contours) and R (blue contours) for a mixture of adiabatic, isobaric 
and isothermal processes, calculated assuming small-amplitude un¬ 
correlated fluctuations of the three processes. The total variance of 
all three processes is normalised to unity. The amplitude of adia¬ 
batic fluctuations increases towards the right, isobaric fluctuations 
upwards and isothermal fluctuations decrease with increasing dis¬ 
tance from the origin. The plot on the top shows the predicted 
values for an unperturbed temperature To = 2 keV and the plot 
on the bottom assumes To = 1.6 keV. These two temperatures 
bracket the range of temperatures observed for M87 in the radial 
range 0.5 — 5'. Boxes schematically show the typical values of R 
and C calculated for several regions in M87. The “arms” and the 
“center” regions are clearly dominated by isobaric fluctuations. On 
the contrary, the circular “ring” is dom inated by adiabatic fluctu¬ 
ations. The outer region (see Sec. 15.31 below! is consistent with a 
mixture of several processes, whith isothermal fluctuations making 
the largest contribution. 


4. TESTING R AND C WITH SIMULATED IMAGES 
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Radius, arcmin 

Fig. 4. — Azimuthally-a veraged radial pr ofile for the deprojected 
temperature in M87 from lChurazov et al.l 12008 *). 

Even in the simplest scenario, when all fluctuations in 
the ICM correspond to only one type of process, other 
complications remain, e.g., radial variations of the mean 
temperature profile, projection effects and, possibly, non¬ 
linear responses to large-amplitude fluctuations. Even 
though the temperature profile is known, the location of 
each individual fluctuation along the line of sight is not. 
The latter effect can be corrected for only in a statis¬ 
tical sense. We proceed by testing the effect of radial 
temperature profiles and different amplitudes of fluctu¬ 
ation on R and C with simulated cluster images using 
the de-projected radia l temperature profile for M87 from 
IChurazov et al.l (120081) . shown in Fig. 0 

The images are produced by generating three indepen¬ 
dent realizations of the density fluctuation patterns (one 
for each of the three processes) as Pi/po = e 5i , where 5 t 
is a Gaussian random field with a given 3D power spec¬ 
trum P, representing fluctuations for process i. This 
log-normal form is used to allow for relatively large am¬ 
plitude density fluctuations. The temperature fluctua¬ 
tions Ti/To are calculated for each process separately, 

e.g. Ti/To = ( pi/po ) 2 ^ 3 for an adiabatic process. These 
values are then used to calculate the fluctuation of the 
volume emissivity ei/eo = The final value of the 

emissivity in each band is calculated as e = eoY^=i 
The resulting 3D cubes are projected onto one plane to 
produce simulated images in each band. For these tests 
we assume, for each process, the same shape of the power 
spectrum P oc k ~, where k is the wavenumber. The 
choice of a k ~ 3 spectrum is rather arbitrary. It has the 
advantage that in 3D it corresponds to the same am¬ 


plitude at all scales and is not far from the canonical 
Kolmogorov slope of fc -11 / 3 , which is c onsistent wit h the 
obser ved spectrum in M87 and Perseus (|Zhuravleva et alJ 

l20lif . 

The soft (1.0-3.5 keV) and hard (3.5-7.5 keV) band 
simulated images are processed in the same way as real 
images (see (0 below), except for the correction for Pois¬ 
son noise. Namely, we fit /3-model to each image, divide 
the image by the /3-model and calculate C and R over a 
range of angular scales. 

In the following two subsections f M.ll and M.2D we per¬ 
form this analysis for pure adiabatic or pure isobaric fluc¬ 
tuations, varying their amplitudes and assuming either 
an isothermal cluster or one with a radial temperature 
profile. Isothermal fluctuations are not considered sepa¬ 
rately, since they have the same impact on the emissivity 
in all bands, so the emissivity patterns and images are 
identical in all bands, except for the overall radial tem¬ 
perature profile which largely cancels out by fitting and 
dividing the images by the projected /3-model. In § 14.31 
we consider a mixture of all three processes. 

4.1. Pure adiabatic fluctuations 

Figure [5] summarizes the results of simulations for pure 
adiabatic fluctuations. Left and right columns corre¬ 
spond to a cluster with a radial temperature profile and 
to an isothermal cluster, respectively. Top and bottom 
rows are for large and small amplitudet|j, respectively. 
The top section in each plot shows the spectra of the 
recovered amplitudes in the soft (red) and hard (blue) 
bands. In the middle sections we plot the correlation co¬ 
efficient C, and the pink line in the bottom panel shows 
the ratio R as defined in eq. [6] The dashed line shows 
the expected value of R for an adiabatic process at an 
unperturbed temperature To = 2 keV. 

For adiabatic fluctuations we can conclude that, for 
these temperatures and energy bands, the correlation co¬ 
efficient C is in the range 0.95-1, and R is between 1.5 
and 2, even for a varying temperature profile and for 
large or small fluctuations. Large fluctuations introduce 
a trend in the value of R, where hard band fluctuations 
increasingly dominate towards smaller spatial scales, but 
their value remains in the range expected for an adiabatic 
process. Therefore, pure adiabatic fluctuations should be 
clearly identifiable in the data by their very high correla¬ 
tion coefficient and large amplitude ratio, C ~ 1, R > 1. 

For a temperature of 2 keV the expected value of R 
for adiabatic fluctuations is ~ 1.6 (see horizontal dashed 
line in Fig. 0. Note that for the temperature profile 
adopted here (see Fig. 0, the temperature is about 1.7 
keV over a significant part of the simulated volume. The 
green curves in the top panel in Fig. [T]show that for this 
temperature, the expected ratio of amplitudes is closer 
to 1.8, in excellent agreement with the R value of these 
simulations. This explains why the recovered value of R 
is larger in the left column of Fig. 0 which uses a more 
realistic temperature profile. 

8 Since our simulated box has finite size, the projected surface 
brightness distribution does not correspond to a perfect /3-model. 
This causes a hump in the amplitude on the largest scales, compa¬ 
rable to the size of the box. 
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Fig. 5.— Test of the recovered C and R values for simulated images with pure adiabatic fluctuations. Left: Average radial temperature 
profile as measured for M87. Right: isothermal cluster at T = 2 keV. The top panels correspond to large amplitude fluctuations and 
the lower panels to lower amplitudes. The amplitude of the fluctuations observed in M87 has intermediate values. In the top section of 
each panel, the blue and red lines show the hard and soft amplitude spectra, respectively. In the middle section, solid blue lines show the 
correlation coefficient C and in the bottom section the solid purple lines show the relative amplitude R. Pure adiabatic fluctuations show 
high correlation coefficient C and high value of R > 1. The dashed lines in the bottom panel mark the expected value of R for adiabatic 
fluctuations around To = 2 keV. 
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4.2. Pure isobaric fluctuations 

The situation is more complex for isobaric fluctuations. 
As can be seen in Fig. [I] the hard band response is close 
to zero in the temperature range of interest, i.e. isobaric 
density fluctuations might show up clearly in the soft 
band but be weak or absent in the hard band. Addition¬ 
ally, the response in the hard band changes from positive 
to negative values as the temperature fluctuates around 
~ 2 keV and this critical temperature also depends on 
the amplitude of the fluctuations. Therefore, the hard 
band response depends sensitively on the amplitude of 
the fluctuations and on the local unperturbed temper¬ 
ature, which can produce correlated, anti-correlated or 
apparently unrelated fluctuations in the soft and hard 
band images. 

Fig. [G] shows the same simulations as in Fig. [5] but for 
pure isobaric fluctuations. In all cases, the hard band 
amplitudes are smaller by a factor of a few than in the 
soft band, and the amplitude ratio is low, as expected. 
Very small fluctuations around Tq = 2 keV (bottom right 
panel) remain in the positive response range and the cor¬ 
relation coefficient remains high and positive. Large am¬ 
plitude fluctuations on the other hand, produce a nega¬ 
tive hard band response, so the fluctuations can be highly 
anti-correlated. For the left panels, which are simulated 
using the radial temperature profile of Fig. [H the ma¬ 
jority of the simulated region is within 6' from the centre 
and therefore has a temperature lower than 2 keV, pro¬ 
ducing mostly anti-correlated hard and soft fluctuations. 

Isobaric fluctuations can be identified more reliably 
through the low amplitude of the hard band fluctua¬ 
tions (|i?| << 1), while the correlation coefficient must 
be interpreted with care, considering the amplitude of 
the fluctuations and the local unperturbed temperature, 
although small (C ~ 0) and negative values of C always 
correspond to isobaric fluctuations. 

4.3. Mixture of three processes 

We now consider a combination of three different types 
of processes by combining the simulated images described 
above and including ‘isothermal” simulated images as 
well, where the fluctuation amplitudes are the same in 
soft and hard bands. To construct the images we simply 
generate three patterns of density fluctuations (modu¬ 
lated by a global /3-model), convert each pattern into 
soft and hard band emissivities, add the flux contribu¬ 
tions from the three processes in each volume element 
and project onto one plane. We measure the correlation 
coefficient C and the ratio R of each pair of images and 
compare these values with the R — C maps, derived for 
an isothermal cluster and linear fluctuations. As before 
we use simulated images in the 1-3.5 keV and 3.5-7.5 
keV bands and vary the relative amplitudes of the three 
processes in independent runs. 

The relative contributions of the three processes in the 
simulated images were chosen to cover several limiting 
cases: each single process dominating, equal contribu¬ 
tions of two processes and equal contributions from the 
three processes. The positions corresponding to these 
setups are marked with crosses in Fig. [7] with single 
processes marked in red, green and cyan. Two equal- 


amplitude processes cases are marked in yellow, black 
and white, and the case where the three processes have 
equal amplitudes is marked in blue. We then processed 
each pair of simulated images as we would do with the 
real data. We calculated R(k re f ) and C(k re f ) values 
for each pair and plotted the resulting values as con¬ 
tours in blue and magenta in Fig. [7] Here we choose 
k re f = 1CT 1 kpc, corresponding to the middle of the 
range of scales we probe in M87. The intersection of 
the contours yields the recovered values of relative am¬ 
plitudes of the adiabatic, isothermal and isobaric pro¬ 
cess. These intersections are marked with ellipses in the 
same color as the corresponding input parameters that 
are marked with crosses. 

The R—C map, that we use to compare the outcome of 
the simulations, was calculated for small-amplitude fluc¬ 
tuations around a constant temperature To = 1.6 keV, 
which is representative of the unperturbed temperature 
range used in the simulations. The results show that, 
even in the presence of a global temperature gradient, 
the leading type of fluctuation is correctly recovered by 
comparing the measured values of C and R with the pre¬ 
dictions for a constant temperature cluster. Indeed, for 
the majority of cases the recovered characteristics of fluc¬ 
tuations (ellipses in Fig. [7]) are close to the input values 
(crosses in Fig. 0. For a combination of two or more 
processes the accuracy deteriorates but the dominating 
processes and their relative fractions can still be qualita¬ 
tively recovered. 

We note that the relative responses of both bands are 
sensitive to the amplitude of fluctuations when the un¬ 
perturbed temperature is low (see low To range in Fig. 
[]]). In this case large amplitude fluctuations can lead to 
R values that fall outside the range of the R — C maps 
in Fig. [3] that was calculated assuming small-amplitude 
fluctuations. Large-amplitude adiabatic fluctuations will 
simply produce larger R values than those predicted 
for the pure adiabatic case (right corner of the R — C 
map), while large isobaric fluctuations can have R values 
smaller than the minimum predicted for small-amplitude 
fluctuations (top left corner of the map). Nevertheless, 
even for large-amplitude fluctuations, adiabatic and iso¬ 
baric fluctuations remain well-separated in the R — C 
map. 

5. NATURE OF THE SURFACE BRIGHTNESS 
FLUCTUATIONS IN M87 

We use archival Chandra observations of M87 with Ob- 
sIDs 2707, 3717, 5826, 5827, 5828, 6186, 7210, 7211, 7212 
and 11783. The d ata were process e d follo wing the proce¬ 
dure described in lYikhlinin et al.l (l2005h . The co-added 
final image has a pixel size of 1" and fairly uniform total 
exposure of about 420 ks within 7' in radius of the cen¬ 
tre, which is the area we will use for the followinganaly- 
sis. The well studi ed X-ray jet fe.g.. lMarshall et alll2002t 
lHarris et al.lI2003H and the point sources, mai nly corre¬ 
sponding to X-ray binaries dJordan et al.ll2004l) , contam¬ 
inate the signal from the ICM fluctuations and were ex¬ 
cised from the images before further analysis. 

We use X-ray images of M87 in the 1-3.5 keV and 3.5- 
7.5 keV bands to calculate the power and cross-power 
spectra of fluctuations, the correlation coefficient C and 












Nature of X-ray Surface Brightness Fluctuations in M87 


9 


k, kpc -1 k, kpc -1 




10- 3 10- 2 10- 1 IO- 3 IQ- 2 10 -i 


k, arcsec -1 


k, arcsec -1 


k, kpc -1 



10 -3 IQ -2 10 -i 


k, kpc -1 



10" 3 IQ -2 10 -i 


k, arcsec 1 


k, arcsec 1 


Fig. 6. — The same as in Fig. [5] but for pure isobaric fluctuations. As expected, pure isobaric fluctuations lead to low values of R < 1 and 
small (or even negative) correlation coefficient C (except for the case of a perfectly isothermal cluster and very small fluctuations). The 
dashed lines in the bottom panel mark the range of values of R expected for isobaric fluctuations around To = 2 keV. 
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adiabatic fluctuation amplitude 

Fig. 7.— Simulated and recovered relative amplitudes of adia¬ 
batic, isobaric and isothermal fluctuations for 7 pairs of simulated 
images. Input values are marked with crosses. Blue and magenta 
lines show the contours, corresponding to C and R value obtained 
from the simulated images. The intersection of these contours de¬ 
fine the recovered values of relative amplitudes of adiabatic, iso¬ 
baric and isothermal fluctuations. These intersections are shown 
with ellipses (same color as the input values). When a single pro¬ 
cess dominates the density fluctuations, the recovered values are 
close to the input values. For a combination of two or more pro¬ 
cesses the accuracy deteriorates but the dominating processes and 
their relative fractions are qualitatively recovered. 

the amplitude ratio R. We fit ^-models to the images in 
each energy band and then divide the images by these 
models to remove the large scale symmetric structure. 
The best fitting /3-model parameters are R c = 0.23' and 
/3 = 0.37 for the 1-3.5 keV image and R c = 0.03' and 
/3 = 0.30 for the 3.5-7.5 keV image. The smaller value of 
/3 for the hard band reflects the global radial temperature 
increase that makes the hard band image less peaked. 
The resulting “flattened” images are shown in Fig. [2] 

We first calculate the amplitude and cross spectrum 
for the entire region, masking out only point sources and 
the jet. This region of the ICM is perturbed by dif¬ 
ferent processes so our aim is to establish whether one 
of these clearly dominates the fluctuations in the clus¬ 
ter. The resulting spectra are shown in Fig. [8] The 
amplitude spectra in the top panel show a very slight 
dependence of amplitude with spatial scale. Perfectly 
horizontal amplitude spectra correspond to a 3D power 
spectrum P 3 d oc k~ 3 so this slight decrease of amplitude 
with k indicates a slightly steeper P 3 d- This is con- 
sis tent with the power spe ctral slope obtained for M87 
by I Zhuravleva et akl (|2014H . although in that paper the 
arms structure was excluded from the analysis making 
the amplitudes smaller , and for the Perseus cluster by 
iZhuravleva et all (12015H . The small R and C values re¬ 
veal the strong contribution of isobaric processes to the 
fluctuation pattern of this region. For spatial scales of 
about 100" we find R = 0.35 and C = 0.5. This value of 
R occurs only at the top left corner of the R — C maps 
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Fig. 8.— Power and cross-power spectra, correlation coefficient 
and ratio R for the 14' x 14' Chandra images of M87 in the 1— 
3.5 keV and 3.5-7.5 keV bands. Top panel: power spectra of 
fluctuation images in the soft (red) and hard (blue) bands. Shaded 
regions show the uncertainty in the measured spectra. The cross¬ 
power spectrum between these bands is shown with the black solid 
line. Middle panel: correlation coefficient C. Bottom panel: 
amplitudes ratio R = P 12 /P\- The dashed lines mark the expected 
value of C for pure isobaric and adiabatic fluctuations, around 
To = 2 keV. The low values of R and C indicate that isobaric 
fluctuations account for 70-80% of the variance in this region. In 
the soft band the total variance is dominated by the contribution 
of filaments/arms of cool gas. 

in Fig. [3l where isobaric fluctuations are responsible for 
about 70-80% of the total variance. In this corner, the 
coherence is very sensitive to the underlying temperature 
so we only remark that, as expected for both maps, the 
absolute value of the coherence should be significantly 
lower than 1 but its precise value is hard to predict. 

We now proceed with the analysis of selected regions 
of the image. We identified structures which might be 
dominated by different processes. From Fig. [2] it is clear 
that the ring structures around 3' from t he centre cor¬ 
res pond to adiabatic fluctuations (see e.g.. lForman et akl 
120071 where this structure is identified as a weak shock), 
while the arms stretching to the east and south-west are 
isobaric, since they do not appear in the hard band im¬ 
age. Guided by these arguments we subdivided the im¬ 
age into three regions: “arms”, “ring” and “outer”. We 
now extract the power spectra of these regions to verify 
if their R — C values match our expectations and estab¬ 
lish the contribution of the different processes in regions 
with less prominent features. 

5.1. Arms 

The filamentary regions extending to about T from the 
centre towards the east and south-west , the arms, were 
identi fied in the X-ray image of M87 bv iBohringer et all 
( 1995 .) who pointed out their spatial correlation with sim¬ 
ilar structures seen in the radio band. In lChurazov et all 
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420011 ) . these structures were explained as gas entrained 
from the core region by rising buoyant bubbles of rela¬ 
tivistic plasma. Later (Molendil (j2002h noted that these 
regions were significantly colder than the surrounding 
medium, and required a two-temperature model for the 
cluster, with tem pera tures as low as 0.8-1 keV for the 
arms. IForman et al.1 (j2005T) produced new temperature 
maps of the M87 cluster confirming the lower tempera¬ 
ture of both arms while finding no strong difference in 
the metallicity of the gas in the arms compared to the 
surrounding medium. 

The part of the M87 image selected for our analysis of 
the arms is shown in the left panel of Fig. [51 The width 
of the selected region is about 130", so we will restrict 
the analysis to wave numbers k > 8 x 10 -3 arcsec -1 . 
The power spectra are shown in the left panel in Fig. 1101 
On the smallest scales, with k > 2 x 10 -2 arcsec -1 , un¬ 
resolved point sources make a significant contribution to 
the variance of the hard band image. Since the spec¬ 
tra of unresolved sources are typically much harder than 
the spectra of the M87 gas, these sources create an up¬ 
turn in the hard band power spectrum and also affect 
the amplitude ratio and correlation coefficient. Since we 
are interested in the arm structures only, we restrict the 
analysis to k < 2 x 10 -2 arcsec -1 . In this range of scales, 
the correlation coefficient and amplitude ratios are small 
with values of R ~ 0.1 and C ~ 0.35, very similar to 
those of the entire image, so the same interpretation ap¬ 
plies here, only shifting the relative importance of the 
isobaric fluctuations to even higher values. The corre¬ 
sponding region in the R — C maps is labeled ‘arms” in 
Fig. [3] 

The value of R ~ 0.2 can be compared with the predic¬ 
tions for isobaric fluctuations plotted in the bottom panel 
in Fig. [1] For an initial temperature of To ~ 2 keV, small 
variations of density (green curve) will produce an am¬ 
plitude ratio of this order of magnitude. Pure isothermal 
(i? ~ 1) and adiabatic processes (i? > 1) can be safely 
excluded, while isobaric fluctuations in these regions are 
completely consistent with the amplitude ratios. 

5.2. Rings at 2!-St.75 (11-17 kpc) 

The ring sections of increased surface brightness at 
radii of 3-3.75' fro m the centre were firs t identified by 
lYoung et~a l. (2002) from Chandra images. IForman et al.1 
(2005) present a temperature map of this region show¬ 
ing an increase in temperature at the rings, compared to 
the immediate surrounding gas, which together with the 
nearly circular appearanc e of the ring su pport its inter ¬ 
pretation as a weak shock iFbr man et al.l (see also [2007). 
Milli on et al.l (120101) used deep Chandra images to con¬ 
firm its weak shock nature. 

The ring appears more prominent in the 3.5-7.5 keV 
band than in the softer bands, consistent with an adi¬ 
abatic origin (top panel in Fig. [[]). We calculated R 
and C for a section of the image covering the northern 
part of the ring, excising only the point sources from 
this region, as shown in the central panel in Fig. [9] The 
thickness of the section including the ring is about 150", 
so we will focus on size scales smaller or equal to this, 
with k > 6 x 10 -3 arcsec -1 . As can be seen in the central 
panel in Fig. [10l around k = 10 -2 arcsec -1 , C > 0.95 and 


R = 1.3. The temperature in this region is about 2 keV, 
so we should interpret the results using the top panel 
in Fig. [3] The location of the intersections of the cor¬ 
responding contours suggests the following distribution 
of amplitudes of the contributing processes: adiabatic ~ 
0.8, isobaric ~ 0.3, isothermal ~ 0.5. In terms of variance 
~ 60% of the total variance can be attributed to adia¬ 
batic processes, ~10% of the variance to isobaric pro¬ 
cesses and ^30% of the variance to isothermal processes. 
For higher mean temperatures, the adiabatic contribu¬ 
tion would increase further. Therefore, adiabatic fluc¬ 
tuations dominate the variance in this region although 
there can be a non-negligible contribution from at least 
one other process. 

5.3. Outer region 

We now proceed with the analysis of the outer region 
shown in the right panel of Fig. [9l where the central 2', 
the arms and the shock have been excluded. The result¬ 
ing power spectra for the 1-3.5 and 3.5- 7.5 keV band 
images are shown in th e right panel of Fig. [10] (see also 
I Zhuravleva et al.ll2014l ). We consider only scales ~100" 
and larger, i.e. below k = 10 -2 arcsec -1 since in this re¬ 
gion the surface brightness is very low in the hard band 
and the power drops an order of magnitude below the 
Poisson noise level slightly above these scales. Around 
k = 10 -2 arcsec -1 , both C and R are about 0.9, point¬ 
ing to a mixture of the three types of processes. For 
mean temperature To ~ 2 keV, the values R ~ 0.9 and 
C ~ 0.9 imply relative amplitudes of ~0.4, 0.4 and 0.8 
produced by adiabatic, isobaric and isothermal processes, 
respectively. In terms of variance, the majority, about 
70%, can be attributed to isothermal processes with the 
remaining ~ 30% distributed equally between adiabatic 
and isobaric processes. For a lower temperature To ~ 1.6 
keV the isothermal processes would have an even higher 
fraction of the total variance, although a temperature of 
1.6 keV is probably too low for this outer region and the 
first estimate is more accurate. Since this region is larger 
than the previous two, we can also focus on fluctuations 
on larger scales. At k = 3 x 10 -3 arcsec -1 , or correspond¬ 
ingly spatial scales of ~ 30 kpc, both the coherence and 
the amplitude ratio drop, pointing to an increase in the 
relative amplitude of isobaric fluctuations. 

Once the prominent structures such as the arms and 
weak shocks are excluded, the outer region appears to be 
dominated by isothermal fluctuations, at least for fluctu¬ 
ations on scales of ~ 10 kpc. In our analysis, isothermal 
structures simply correspond to density fluctuations that 
produce equal amplitude fluctuations in the soft and hard 
band fluxes, such as, e.g., X-ray cavities produced by ra¬ 
dio bubbles, or the global asymmetry of the cluster. The 
dominance of isothermal fluctuations of these scales in 
this region is c onsistent with the presence of several b ub- 
bles (see, e.g., lOwen et~Mll2000HForman et al.ll2007t ). In 
addition to these bubbles, at distances of about 6'.7 from 
the centre the gas is likely disturbed by large scale slosh¬ 
ing ( e.g. Fig. 2 in Sim ionescu et al .1 [20 1 (1 [Werner et al.1 
l2015f) .The most prominent cold front presumably asso¬ 
ciated with sloshing is located ~20' north of M87. This 
feature is outside the region studied here. There is also 
a feature ~6' south of M87 (green area in the lower left 
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Fig. 9.— Selected image regions. Left: arms, center: shock, right: rest. The sections correspond to the same 1—3.5 keV image of M87, 
divided by the best-fitting /5-model. All images are on the same scale. The size of the largest image is 14'x 14'. 
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Fig. 10.— The same as in Fig. [8] for three separate regions shown in Fig. [9] Left panel: the region dominated by cool arms. The 
value of R is close to that expected for isobaric fluctuations around mean temperature To = 2 keV. The correlation coefficient is small, as 
expected. Central panel: the shock region. Two dashed black lines correspond to pure adiabatic and isobaric processes for 3 keV gas. 
The values of R and C suggest that adiabatic process (weak shock) dominates. Rights panel: ”outer” region (no arms, no central region, 
no shock) is only robust for k < 10“ 2 arcsec —1 since there are very few counts in this region and the hard band power drops by an order of 
magnitude below the Poisson noise level. Around and below k = 10 -2 arcsec -1 both the correlation coefficient and ratio are below unity 
(see text for details). 


quadrant of the right panel in Fig. 0 that might be a 
fainter counterpart of the more prominent northern cold 
front, similar to the structures predicted by 2D numerical 
simulations of gas sloshing, s hown in the bottom panel 
in Fig. 11 in lChurazov et al] (120031 ) . 

5.4. Central 2! 

The central 2' (radius) area appears disturbed by sev¬ 
eral processes so we study it separately. Fig. [Tl] shows 
this region in three energy bands: 0.5-1, 1 3.5 and 3.5- 
7.5 keV. Each image has been divided by its best-fitting 
/3-model to remove the symmetric large scale structure 
and highlight the fluctuations, which are shown on the 
same colour scale. Almost all the bright features, espe¬ 
cially the base of the arms toward the east and south¬ 
west, appear stronger in the softest band, indicative of 
their isobaric nature. A few areas, marked by the red 
arrows, however, appear brighter in the 1 -3.5 band than 
in the 0.5-1 keV band, which points towards adiabatic 


fluctuations. These regions include the “bud” structure 
identified bv lForman et al.1 (|2005l [2007h as a budding ra¬ 
dio bubble. In the hard 3.5-7.5 keV band image (right 
panel of Fig. m, only this inner ring structure is visible, 
suggesting an adiabatic nature of the rings and corrobo¬ 
rating the isobaric nature of the base of the arms. 

The central region is cooler than the rest of the clus¬ 
ter, as can be seen in the temperature profile in Fig. [2 
averaging about 1.6 keV. For this low temperature, there 
are noticeable differences in the response to density fluc¬ 
tuations in the 0.5-1 keV and 1-3.5 keV bands, which 
we can exploit to study the origin of the fluctuations 
further. Figure [T3] shows the amplitude ratios for adia¬ 
batic and isobaric density fluctuations, similarly to Fig. 
[Qbut this time comparing the 0.5-1 keV and 1-3.5 keV 
bands. For unperturbed temperatures below 2 keV, the 
fractional amplitudes deviate strongly from unity, espe¬ 
cially for isobaric processes. 
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Fig. 11. — Central, 1' radius region, from left to right: 0.5-1. 1-3.5 and 3.5-7.5 keV images, each divided by its best-fitting /3-model. 
Note how the fluctuations corresponding to the base of the arms decrease with increasing energy while the ring structure (marked with red 
arrows in the central plot) becomes more dominant with respect to other feature with increasing energy. These trends are characteristic of 
isobaric and adiabatic fluctuations, respectively. 
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Fig. 12. — Amplitude ratios for 0.5—1 keV and 1—3.5 keV band 
images for adiabatic (top panel) and isobaric (bottom panel) pro¬ 
cesses, respectively. The fluctuations are modelled as density vari¬ 
ations from an initial value p 0 by a factor G (p = p 0 * G). At 
temperatures below 3 keV, these bands have amplitude ratios sig¬ 
nificantly different from unity. The adiabatic fluctuations appear 
significantly larger in the hard band, while isobaric fluctuations 
appear larger in the soft band. The thick green lines in each plot 
correspond to small amplitude fluctuations, while the thin lines 
correspond to larger (nonlinear) fluctuations as labeled in the fig¬ 
ure. 


The power spectrum analysis of the central region 
for the very soft 0.5-1 keV and soft 1-3.5 keV energy 
bands is shown in Fig. [13] The correlation coefficient 
is C ~ 0.9 — 0.7 and i? ~ 0.5 — 0.4. These value can 
be compared directly with the curves in Fig. 1121 For 
an unperturbed temperature To = 2 keV, an amplitude 
ratio of 0.5 requires a factor of ^2 density increase due 
to isobaric processes. If the unperturbed temperature 
is Tq = 1.6, then the required density increase ~ 1.5 is 
more modest. The correlation coefficient drops slightly 
towards smaller spatial scales (higher values of k) hinting 
at uncorrelated fluctuations of a different type becoming 
stronger at smaller scales, probably the adiabatic inner 
ring structure marked with red arrows in Fig. 1111 
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Fig. 13. — Top: Amplitude (red, blue) and cross spectrum (black) 
between 0.5-1 keV and 1-3. 5 ke V images, for the central region (left 
and central images in Fig. ED. Middle: correlation coefficient, C. 
Bottom: amplitude ratio R = Pi^/Pi- R = 0.5 is consistent 
with mainly isobaric fluctuations, where the density is lower than 
ambient by 10% if the ambient temperature is To = 1.2 keV (green 
line in Fig. I12D , by a factor of 1.5 from an ambient temperature 
of 1.8 keV or by a factor of 2 from a ambient temper atur e of 2 
keV (green, blue and black lines, respectively, in Fig. I12D . The 
correlation coefficient C ranging from 0.7 to 0.9 is consistent with 
dominant isobaric fluctuations, contaminated by other processes. 

For this region, we have also calculated the values of R 
and C for the images in the 1 3.5 and 3.5-7.5 keV energy 
bands (see Fig. nut. On scales ~2' (or equivalently ~9 
kpc), similar to the radius of this central region, C ~ 0.4 
and R ~ 0.3, which correspond to the top left corner 
of the R — C maps in Fig. [3] i.e. mostly isobaric fluc¬ 
tuations, confirming the result of the low energy bands. 
On scales of ^20", about the thickness of the ring in the 
right panel in Fig. [TIJ the values of R and C both in¬ 
crease to ~ 0.5, corresponding to an increased contribu- 
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Fig. 14.— Top: Amplitude (red, blue) and cross spectrum (black) 
between 1—3.5 keV and 3.5—7.5 ke V im ages, for the central region 
(central and right images in Fig. 1 11 1) . Middle: correlation coef¬ 
ficient C. Bottom: amplitude ratio R = Pi^/Pi- On scales of 
2'or equivalently 9 kpc, C = 0.4 and R = 0.3, which correspond 
to the top left corner of the R — C maps in Fig. [3] i.e. mostly 
isobaric fluctuations. On scales of 2 0", about the thickness of the 
ring in the right panel in Fig. 1111 the values of R and C both 
increase to 0.5, corresponding to an increased contribution from 
adiabatic processes which make up about 30% of the variance of 
density fluctuations at these smaller scales. 

tion from adiabatic processes which make up about 30% 
of the variance of density fluctuations at these smaller 
scales. 

6. ENERGY CONTENT OF THE FLUCTUATIONS 

We now relate the total excess energy, associated with 
fluctuations, to the observed integrated variance M 2 . 
Consider, for example, the density fluctuations produced 
by a collection of bubbles of relativistic plasma, occu¬ 
pying a fraction /j < 1 of the cluster volume. We as¬ 
sume that the bubbles are completely devoid from ther¬ 
mal plasma. Therefore, in X-ray observations they corre¬ 
spond to density fluctuations Sp/(p) ~ —1. The resulting 
total variance (to leading order in fb) is 



Since the free energy, as sociated with bubble s is their 
enthalpy Eb = ^zjPVb ( Churazov et al . 200 111 , we can 
write Eb = 7/bAthermab where Vb is the volume occupied 
by bubbles and ^thermal = ■^ijPV is the total energy in 
the volume of the region under consideration, V. Thus 

T7 — -7 M 2 . (8) 

-^thermal 

Similar expressions can be derived for density fluctua¬ 
tions due to sound waves or gravity waves, they all have 


the form • E g. uctuation ~ «M 2 , where a is of order unity 
(Zhuravleva et al. in prep). Therefore, independently of 
the nature of fluctuations, we can crudely estimate the 
total energy associated with fluctuations as « M 2 . 

To estimate the total variance of density fluctuations, 
M 2 , we will use the power or amplitude spectra of the 
1-3.5 keV images, since these correspond most closely to 
density fluctuations. In terms of the power this reads 

M 2 = [ 2 P 3D (k)4Trk 2 dk (9) 

Jfci 

where P 3 o(k) is the de-projected power spectrum and k 
is the spatial wave number. In terms of the 3D amplitude 
of fluctuations, as plotted in the top panels in Figs. [8]and 
EH this equates to 

M 2 = [" 2 ^PP-dk. (10) 

Jk ! k 

We integrated the 3D amplitude spectra obtained for 
the entire 14' x 14' region in the 1-3.5 keV band, between 
k\ = 0.005 arcsec -1 and fc 2 = 0.05 arcsec -1 obtaining a 
value of M 2 ~ 0.05. This limited range in spatial scales 
used for the integration (see Eq. m evidently results in a 
lower limit to the total variance of density fluctuations. 
We do not expect a large contribution to the variance 
from larger k since the amplitude spectrum declines with 
increasing k. On the other end, k < fci, corresponding 
to larger scales, could potentially provide a significant 
contribution to the variance. However, the value of k\ is 
not far from 1/R where R is the size of the region. It is 
therefore unlikely that our estimates of the total variance 
at scales < R are underestimated by a large factor. 

For the 14' x 14' region, M 2 ~ 0.05 as discussed above, 
implying that in this region the amount of non-thermal 
energy is ~ 5% of the total thermal energy. From here 
we can estimate that the time fdiss needed to convert 
this non-thermal energy into heat in order to balance the 

Cooling losses is tdiss = ^cooling x -Efluctuation/E/thermal ~ 

0.05 tcooiing- For typical conditions in cool-core clusters, 
this corresponds to the dynamic timescale, within a fac¬ 
tor of a few. Th is timescale is consistent wi th the AGN 
feedback model (|Churazov et al.l 120001 . l200Hl when bub¬ 
bles of relativistic plasma initially capture much of the 
AGN mechanical power and then transfer energy to the 
gas on times scales set by the rise of buoyant bubbles. 
According to this model, the energy goes into turbulence 
and gravity waves —including the contribution of en¬ 
trained gas— and dissipates into heat. 

A similar analysis can be done for the amplitude spec¬ 
tra of selected regions, as shown in Figs. liilandfTOl In this 
case we took two approaches. First, we integrate the am¬ 
plitude spectrum directly, in a range of scales where this 
spectrum is reliably measured. Then we used this range 
of k to fit a powerlaw model to the amplitude spectra and 
integrated the model over the same range of k as used 
above for the entire image i.e. k = 0.005 — 0.05arcsec -1 . 
The results, with both methods in each region, are sum¬ 
marised in Table E The M 2 in the arms region is about 
four times higher than for the entire image, which has 
the arms contribution to the power diluted by the less 
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Nature of X-ray Surface Brightness Fluctuations in M87 


Region 

k range 


k range 

M m 

all 

0.005-0.05 

0.058 

0.005-0.05 

0.053 

center 

0.007-0.1 

0.11 

0.005-0.05 

0.098 

arms 

0.01-0.1 

0.18 

0.005-0.05 

0.22 

shock 

0.007-0.04 

0.016 

0.005-0.05 

0.019 

outer 

0.005-0.02 

0.012 

0.005-0.05 

0.015 


TABLE 1 

Integrated power of density fluctuations for the regions noted in 
the left column. Direct integration of the amplitude spectra are 
quoted in the third column and analytic integration of the 
powerlaw fit to the spectra are quoted in the fifth column. The 

frequency ranges used for each integration are quoted in the 
second and fourth columns. 

perturbed part of the image. The shock and the outer 
regions have much lower amplitudes, with M 2 about an 
order of magnitude lower than the arms region. 

We can estimate the contribution of each section of the 
image to the total power as follows. The outer region, 
where the shocks and arms are excluded, is probably rep¬ 
resentative of the average conditions in the cluster, i.e. 
similar bubbles and other structures should also exist in 
the image regions that have been masked out, perhaps 
at different position along the line of sight than the arms 
and shocks. Therefore, the contribution to the total M 2 
of these mainly isothermal structures is likely similar to 
the M 2 ~ 0.015 contribution found in the outer region. 
The shock region is probably representative of the entire 
ring with inner radius 1'.3 and outer radius 4'.5 from 
the centre, not just the section used for our calculations, 
that masked out the overlying arms and other structures. 
Therefore, if no other shocks are present, the shock will 
contribute with about 0.02 x0.3 = 0.006 to the total M 2 
where the factor of 0.3 accounts for the ratio of areas of 
the shock ring to the entire 14 x 14 image. For the arms 
region, if these are the only such structures in the to¬ 
tal image then their contribution only needs to be scaled 
down by the ratio of areas so that the contribution of 
arms to the total M 2 is ~ 0.22 x 1/9 ~ 0.024. Finally, 
the central region, scaled down by the area ratio con¬ 
tributes with ~ 0.1 x 1/15 = 0.006, which in total add 
up approximately to the variance of the entire region, i.e. 
M 2 ~ 0.05. 

Note that the dissipation timescale discussed above 
corresponds to an average time over the entire 14 x 14' 
region. Isolating the central region, within 2' from the 
centre, the corresponding M 2 ~ 0.1 so here fdiss ~ 
0-1 tcooiing- Repeating these calculations for radii up to 
4' also produces fdiss ~ 0.1 f C ooiing- 

7. CONCLUSIONS 

In this paper, we demonstrate how adiabatic, isobaric 
and isothermal fluctuations in the ICM have different 
impacts on cluster images in different X-ray bands. For 
our chosen energy bands (1-3.5 and 3.5-7.5 keV) adia¬ 
batic fluctuations, such as sound waves or weak shocks, 
produce larger amplitude in the hard band. On the con¬ 
trary, isobaric fluctuations, characteristic of subsonic gas 
displacement, produce larger surface brightness fluctua¬ 
tions in the soft band, while the hard band image might 
show only a small or even a negative response to the 
same fluctuations. Isothermal fluctuations, such as bub¬ 
bles and large scale asymmetries, appear equally strong 


in both bands. We have exploited this behaviour of the 
cluster X-ray emission to study the origin of the features 
in the ICM that produce the observed surface brightness 
fluctuations. Of course, there are a number of caveats, 
associated with this analysis. These include the uncer¬ 
tainty in removal of the global cluster profile using a 
simplistic /3-model; the assumption of small-amplitude 
linear fluctuations, and the assumption that processes of 
different type are spatially uncorrelated. We neverthe¬ 
less believe that the analysis captures important obser¬ 
vational signatures of different processes responsible for 
the observed fluctuations and allows us to identify the 
most prominent ones. 

The most prominent features of the M87 X-ray im¬ 
ages had been identified as cool arms and quasi-sp herical 
shocks (see figs. 2a, 9a, 11a in lForman et a 1.11200 7;. We 
isolated the regions that contain these structures and cal¬ 
culated power and cross-power spectra, the correlation 
coefficient C and the ratio of amplitudes R. The derived 
values of R and C show an excellent agreement with the 
predictions of isobaric fluctuations in the arms region and 
adiabatic fluctuations in the shock region. 

The rest of the cluster is not clearly dominated by a 
single type of fluctuation. We divided the remaining area 
into a central 2' circle and an outer region. In the centre, 
the leading type of fluctuation is isobaric, dominating the 
variance on scales of 100". At smaller scales of 20", adi¬ 
abatic fluctuations appear and account for about 30% of 
the variance. The dominance of isobaric structures in the 
central region is further confirmed by the comparison of 
the soft 1 3.5 keV to the very soft 0.5-1 keV images (see 
left and centre panels in Fig. ED- The amplitude ratio 
R between these bands is about 0.5 which is consistent 
with a factor of ^2 overdensity region in pressure equi¬ 
librium with ambient gas at about 1.6 keV, as shown by 
the predicted amplitude ratios for isobaric fluctuations 
between these bands in Fig. [12] 

For spatial scales of 5-10 kpc, the outer region is dom¬ 
inated by isothermal fluctuations, which include bubbles 
and deviations from the splrerically-symmetric model 
used to remove the large scale gradient produced by the 
potential of the cluster. The measured R and C val¬ 
ues correspond to about 70% of the variance produced 
by isothermal processes and the remainder is distributed 
about equally between adiabatic and isobaric processes. 
Fluctuations on larger scales, of about 30 kpc point to 
a larger contribution from isobaric fluctuations. Part of 
these isobaric and apparently adiabatic fluctuations can 
be unrelated to the central AGN activity, for example as 
a result of gas sloshing produced by in-falling substruc¬ 
ture. 

We show that the energy content of the density fluc¬ 
tuations is proportional to their integrated variance, so 
that by integrating the amplitude spectra of the 1-3.5 
keV images, we can estimate the ratio of non-thermal 
(fluctuations) energy to the thermal energy. Using 
this argument ,we find a non-thermal energy fraction 
-Fluctuation/-^thermal ^ 0.05 foi the entire 14 x 14 re¬ 
gion centered on M87. From this value we estimate the 
dissipation time of the fluctuations as fdiss ~ 0.05f coo iing 
in order for these fluctuations to balance radiative cool¬ 
ing. The total variance depends on the region of the 
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image considered, in particular, in the central 2' radius 
and the central 4' radius regions the integrated variance 
is M 2 ~ 10%, indicating shorter dissipation timescales 
of tdiss ~ O.ltcooiing. These results broadly agree with an 
AGN feedback model mediated by bubbles of relativistic 
plasma. 
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